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ABSTRACT 

We present three Far Ultraviolet Spectroscopic Explorer (FUSE) observations of the Narrow-Line Seyfert 1 
galaxy NGC 405 1 . The most prominent features in the far-ultraviolet (FUV) spectrum are the O VI emission 
and absorption lines and the H I Lyman series absorption lines which are detected up to the Lyman edge. We 
also identify weak emission from N III, C III, and He TL The C III line shows absorption while none is detected 
in the N III and He II lines. In H I and C III we detect two main absorption systems at outflow velocities 
of -50 ±30 and -240 ±40 kms" 1 , as well as a possible third one at ~ —450 km s 1 . These systems are 
consistent in velocity with the 10 absorption systems found previously in C IV, N V, and Si IV, though the 
individual systems are blended together in the FUV spectrum. We estimate column densities of the two main 
absorption systems and find that the H I column density is lower for systems with larger outflow velocity. 
We detect no flux or spectral variations of NGC 405 1 at FUV wavelengths during three epochs spanning one 
year. This is consistent with the optical light curve which shows no variations between the three epochs. It is 
also consistent with the X-ray light curve which shows consistent flux levels at the three epochs of the FUSE 
observations, although the X-ray light curve shows strong variations on much shorter timescales. 
Subject headings: galaxies: active — galaxies: individual (NGC 4051) — galaxies: nuclei — galaxies: Seyfert 
— techniques: spectroscopic — ultraviolet: galaxies 



1. INTRODUCTION 

The Narrow -Line Seyfert 1 (NLS1) galaxy NGC 4051 is 
one of the best-studied active galactic nuclei (AGN) across the 
electromagnetic spectrum due to its brightness and proximity 
(V w 13.5; z = 0.002295 ± 0.000043 from optical emission 
lines; de Vaucouleurs et al. 1991). It has been part of many 
studies since it was first mentioned by Hubble (1926) and 
classified by Seyfert (1943). Collinge et al. (2001) recently 
presented the first high-resolution X-ray and ultraviolet (UV) 
spectra of NGC 405 1 . The Chandra High-Energy Transmis- 
sion Grating Spectrometer (HETGS) spectrum reveals ab- 
sorption and emission lines from hydrogen-like and helium- 
like ions of O, Ne, Mg, and Si. Two distinct blueshifted 
absorption systems were detected: a high-velocity system at 
-2340 ± 130km s" 1 and a low- velocity system at -600 ± 130 
kms~'. In a Chandra Low-Energy Transmission Grating 
Spectrometer (LETGS) observation taken ~ 2 yr later, van 
der Meer et al. (2003) claim to detect a new X-ray absorption 
system at ~ —4500 km s while the —2340 km s _1 system is 
barely visible. An even higher velocity outflowing absorption 
system at ~ —6500 krns" 1 or ~ —16500 kms -1 (depend- 
ing on the identification of an absorption line at 7.1 keV in 
an XMM-Newton EPIC observation) was suggested recently 
by Pounds et al. (2003). The UV spectrum taken with the 
Space Telescope Imaging Spectrograph (STIS) shows strong 
absorption mainly in the lines of C IV, N V, and Si IV. Ten 
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intrinsic UV absorption systems are seen with velocities be- 
tween —650 and 30 km s~ 1 and FWHMs ranging from ~ 30 to 
~ 160 kms -1 . While the low-velocity X-ray absorption sys- 
tem is consistent in velocity with many of the UV absorption 
systems, the high-velocity X-ray absorption seems to have no 
UV counterpart. 

NGC 4051 shows rapid and large-amplitude X-ray flux 
variability by up to an order of magnitude on timescales of 
hours (e.g., Marshall et al. 1983; Uttley et al. 2000; Collinge 
et al. 2001; McHardy et al. 2003; and references therein). 
Uttley et al. (2000) found a strong correlation between the 
variability in the extreme UV (124-188 eV) and X-ray (2- 
10 keV) bands, and they suggested this to indicate that both 
bands are sampling the same power-law continuum. Shem- 
mer et al. (2003) find the variable optical continuum flux to 
correlate with the variable X-ray flux. The optical continuum 
leads the X-ray by 2.4 ± 1.0 days. They also found that the 
variation amplitude in the optical is much smaller than in the 
X-ray, as seen in other NLS Is. 

In this paper we bridge between the X-ray and near UV 
high-resolution spectra by presenting the first Far Ultraviolet 
Spectroscopic Explorer (FUSE) spectra of NGC 4051. The 
FUSE spectra cover the 900-1 180 A band with a resolution 
of ~ 0.06 A (corresponding to ~ 20 km s _1 ) . Our main pro- 
posed goals were ( 1 ) to study the FUV kinematic counterparts 
to the known X-ray and UV absorption lines, (2) to study any 
FUV absorption-line variability, and (3) to relate any FUV 
variability to that seen at X-ray and other wavelengths. The 
observations and data reduction are detailed in § 2. In §§ 3 
and 4 we present the temporal and spectral analyses, and in 
§ 5 we discuss our results. Finally, §6 presents a summary of 
our main results. 

2. OBSERVATIONS AND DATA REDUCTION 

NGC 4051 was observed with FUSE (Moos et al. 2000; 
Sahnow et al. 2000) during three epochs, spanning one year, 
which are listed in Table [l] The observations were carried 
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TABLE 1 

FUSE Observation Log for NGC 405 1 



Obs. ID 


Date and UT 


N m bits d 


Time (ks) b 


B0620201000 


2002 Mai' 29 01:58 


21 


28.7 


C0190101000 


2003 Jan 18 11:43 


16 


34.7 


CO 190 102000 


2003 Mai' 19 07:46 


17 


28.8 



"Number of consecutive orbits. 
b Total effective exposure time. 



out using the 30" x 30" low-resolution (LWRS) aperture in 
the standard observing mode. Each observation consists of 
several consecutive orbits, and in each orbit data from eight 
detectors (each covering a different wavelength range) were 
accumulated. We reduced the raw data using the most recent 
FUSE software (CalFUSE v2.4.0 and the FUSE IDL tools 
version of 2002 July). For each observation we checked the 
data from different detectors and orbits for consistency and 
found no significant deviations in flux. Hence, we used the 
FUSE_REGISTER tool to combine the data sets from the dif- 
ferent detectors and orbits into one spectrum per observation. 

In the following analysis we use the line templates from 
Feldman et al. (2001) to identify the airglow geocoronal emis- 
sion lines and Sembach (1999) to identify interstellar Galactic 
absorption lines. 

3. TEMPORAL ANALYSIS 

Visual comparison of the spectra from the three epochs sug- 
gests striking similarity between them with no apparent time 
variation of the intrinsic spectrum. To assess any time varia- 
tions quantitatively, we computed (/i — fa ) / J (of + ) f° r 
each pair of spectra, where f\ and fa are the flux densities 
at a given wavelength and G\ and Oi are their uncertainties. 
Plots of this quantity versus wavelength revealed that the only 
variable regions in the spectra (including regions where ab- 
sorption lines are present) could be explained by variable geo- 
coronal emission lines or uncertainties in calibration. We also 
measured the average flux densities of the three spectra in 
several wavelength bands which are free from emission and 
absorption; these measurements are listed in Tableland con- 
firm that there are no variations within the ~10% uncertainty 
of each measurement. A careful look at Table|2]suggests that 
the flux of the 2002 Mar 29 spectrum may be somewhat lower 
(by ~ 3%) than for the two spectra taken about a year later, 
especially longward of ~ 990 A. However, this possible vari- 
ation is not highly significant (< la). Since we found the 
spectra from the three epochs to be consistent, we combined 
them into one 92.2 ks spectrum to increase the S/N of our 
analysis (see Figure^. 

We also checked for flux variations within each observa- 
tion by comparing the spectra between single orbits and by 
combining 4-5 orbits together. The spectra from single orbits 
have very poor S/N (of order 1), and no variability is detected 
between them within the limiting S/N of the data. Analyses of 
the spectra made by combining 4-5 orbits together show that 
any FUV variability is less than ~ 20% (considering several 
wavelength bands as given in Table 2). 

We compared the FUSE observations to the long-term 
RXTE 2-10 keV light curve of NGC 405 1 (Lamer et al. 2003; 
McHardy et al. 2003; P. Uttley and I. M. McHardy 2003, 
private communication). During the year of the FUSE obser- 



TABLE 2 

Average Flux-Density Measurements 



Wavelength range 


2002 Mar 29 


2003 Jan 18 


2003 Mar 19 


926.5-927.2 


13.39± 1.22 


13.86± 1.21 


13.37 ± 1.23 


939.8-941.8 


11.01 ±0.76 


10.95 ±0.70 


10.29±0.73 


956.9-958.0 


13.09± 1.09 


12.61 ±0.98 


13.18 ± 1.07 


967.1-968.2 


9.28 ± 1.10 


10.69 ± 1.02 


8.92 ± 1.05 


998.2-999.2 


13.93 ±1.04 


16.81 ±0.94 


15.24±1.01 


1019.2-1020.7 


15.45 ±0.72 


16.92 ±0.64 


16.56 ±0.69 


1030. 1-103 1.6 C 


7.81 ±0.62 


8.31 ±0.56 


8.15±0.59 


1042.0-1043.0 d 


21.38±0.90 


21.69±0.79 


21.33±0.87 


1056.1-1057.6 


12.76 ±0.76 


14.36 ±0.67 


13.16±0.71 


1069.0-1070.0 


13.29 ±0.97 


14.42 ±0.85 


14.50±0.92 


1089.0-1090.0 


12.22 ±1.65 


13.33 ±1.33 


11.29± 1.34 


1095.0-1096.5 


9.98 ±0.87 


11.88±0.72 


11.69±0.75 


1114.4-1116.4 


10.79 ±0.37 


11.97±0.33 


10.99 ±0.35 


1127.4-1129.4 


9.59 ±0.36 


11.41 ±0.33 


10.75 ±0.35 


1145.4-1147.4 


8.97 ±0.41 


9.17±0.34 


9.50±0.38 


1181.3-1183.3 


13.30±0.96 


14.04 ±0.80 


13.56±0.88 


"In units of A in the rest frame. 



b In units of 10 15 erg cm 2 s 'A 
Measurement in the trough of O VI A 1032. 
d Measurement on the red wing of O VI A 1038. 



vations, NGC 4051 was observed with RXTE every two days 
allowing fairly contemporaneous comparison of the X-ray and 
FUV fluxes (although, as noted in § 1, NGC 4051 is highly X- 
ray variable on even shorter timescales). The fluxes from the 
three RXTE observations closest in time to the FUSE observa- 
tions are consistent [the 2-10 keV fluxes are measured to lie 
in the range (3-3. 5) xl0~ 11 ergcm _2 s -1 ]. We also compared 
the FUSE observations to optical observations taken on simi- 
lar dates (within about 10 days) at the Wise Observatory. The 
optical observations show no variations either. This is consis- 
tent with the results of Peterson et al. (2000) and Shemmer 
et al. (2003) who find NGC 4051 to vary by only a few per- 
cent in the optical. 

4. SPECTRAL ANALYSIS 

The most prominent intrinsic features seen in the FUSE 
spectrum (Figure Q are the O VI emission and absorption 
lines and the H I Lyman series absorption lines, which are de- 
tected up to the Lyman edge. We also identify weak intrinsic 
emission from N III, C III, and He II. The C III line shows ab- 
sorption, while no intrinsic absorption is detected in the N III 
and He II lines. All absorption lines are blueshifted relative 
to the systemic velocity, and some show multiple absorption 
systems. 

To study the intrinsic O VI absorption, we fitted each of 
the two doublet emission lines with two Gaussians, keeping 
the width of each Gaussian the same for both doublet lines 
and keeping the flux ratio the same as the oscillator-strength 
ratio (2:1). The fit results are shown in Figure [2] and listed 
in Table [5] The observed spectrum was then divided by this 
emission model, and the resulting normalized spectrum was 
used to create the absorption velocity spectra shown in Fig- 
ure |3jt. The absorption-line profiles of the two O VI lines are 
similar with no noticeable features in the trough which ranges 
from about 50 kms -1 to —800 kms~'. The overall spread 
of velocities is in agreement with the 10 absorption systems 
identified in the UV lines of C IV and N V (see Figure [3j>); 
however, no individual absorption systems are seen in O VI. 
Since the resolution of the FUV spectrum is high enough to 
resolve these systems at least in part (see Figure[3j;), and since 
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FIG. 1 .— The combined 92.2 ks FUSE spectrum of NGC4051 binned to ~ 0.1 A pix . Intrinsic emission lines are marked at their theoretically expected 
positions in red above the spectrum. Airglow geocoronal emission lines and interstellar Galactic absorption lines are marked with ©. The intrinsic H I absorption 
lines are marked at their theoretically expected positions up to the ionization energy in the top panel. 



the troughs of the two lines have similar depths (thus depart- 
ing from the oscillator-strength ratio), we conclude that the 
O VI absorption lines are saturated and have broad enough 
damping wings to create a single trough. The O VI trough is 
not completely black which indicates partial covering of this 
absorber. 7 From Figure [3^ we estimate the absorber covering 
factor to be -0.9-0.95. 

In addition to its FUV lines, the ion O VI is also expected 
to have X-ray absorption lines at 19.135, 19.341, 21.79, and 
22.01 A (E. Behar 2003, private communication; see also 
Pradhan 2000). The O VI AA1032, 1038 absorption sys- 
tems detected in the FUSE spectrum have a combined width 
of — 700 kms -1 . The velocity resolution of the Chandra 
HETGS (LETGS) spectrum, at -20 A, is -350 (600) km s" 1 ; 
hence the O VI absorption lines can be resolved by these in- 
struments. Examination of the HETGS and LETGS spec- 

7 We rule out the possibility that the trough is not black due to an instru- 
mental scattering problem since the Galactic absorption line at ~ 1034 A 
reaches zero flux. We also rule out the possibility of scattering of nuclear 
light into the line of sight sinc e other absorption lines reach zero flux; see the 
discussion at the end of § 15.11 




1010 1020 1030 1040 
Rest Wavelength [A] 



1050 



FIG. 2. — Fit to the O VI emission lines using two Gaussians for each line 
(dashed curve) plotted over the FUSE spectrum (solid curve) which is binned 
to ~ 0.1 A pix -1 . The lo~ uncertainties for the points are also shown. The 
fit parameters are given in Tablel3l 
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TABLE 3 
O VI FIT PARAMETERS 



Gaussian 


Normalization" 


o" b [A] 


FWHM b [kms" 1 ] 


I 


0.540 


1.074 


735 


n 


0.179 


5.831 


3970 



NOTE. — All Gaussian centroids were set to the same velocity. In the rest- 
frame spectrum of the source they correspond to 1031.22 A and 1036.91 A, 
i.e., the fitted O VI emission lines are blueshifted by —200 km s~ 1 relative to 
the optical lines. 



a The O VI A 1032 line normalization in units of 10 13 erg cm 2 s 1 A '. 
The O VI A 1038 line normalization was set to half the O VI A 1032 line 
normalization. 

b o~ and FWHM are the same for the two lines of the doublet. 



tra shows that the data are consistent with these lines being 
present, though the poor S/N prevents any significant detec- 
tion. 

In Figure|5|; we compare the absorption profile of C III with 
that of C IV A 1550 from Collinge et al. (2001). C III shows 
absorption in the velocity ranges of —150 to 50, —330 to 
—230, and possibly around —450 kms . Within the limita- 
tions of the FUV spectral resolution (~ 20 kms -1 ), the main 
absorption systems that were identified in the C IV STIS pro- 
file are also seen in the C III FUSE profile. The recession 
velocity of NGC4051 is ~ +690 kms -1 . Any absorption at 
an apparent rest-frame blueshifted velocity of ~ —690 km s~ 1 
is probably Galactic or contaminated by Galactic absorption, 
and no useful information about the intrinsic absorption in 
NGC 405 1 can be deduced from absorption near this velocity. 

Velocity spectra for the first eight lines in the H I Lyman 
series are shown in Figure |4] The Lya absorption observed 
with STIS (from Collinge et al. 2001) is shown in Figure^ 
together with Lyj3 from FUSE. Lya does not show individ- 
ual intrinsic systems but rather has a broad absorption trough 
which reaches up to a velocity of ~ —1200 kms -1 . All the 
Lyman series lines observed with FUSE show geocoronal Ly- 
man emission lines at ~ —550 to —750 km s . The absorp- 
tion trough of Lyj3 is further contaminated with geocoronal 
O I lines which mask the absorption in the range ~ 90 to —270 
km s~ 1 . The Lyy and Ly 8 lines are shown in Figure |4p. Both 
lines are contaminated at velocities < —430 kms . 

The H I Lyman absorption lines are seen all the way to 
the Lyman edge (see Figure ^ top panel and Figure |5}- We 
identify two main H I Lyman absorption systems: one at 
—50 ± 30 kms -1 which is detected in the lines from Lyy to 
Ly^ (after which the lines are blended together and blended 
with airglow emission lines; see Figure[5}, and one at —240 ± 
40 kms -1 which is detected in the lines from Lyy to LyX 
(after which it is weak and blended with the —50 kms -1 ab- 
sorption lines). All lines identified in the —50 kms -1 system 
are saturated. 

We used the XVOIGT software (Mar & Bailey 1995) to es- 
timate the column densities of the two H I Lyman absorption 
systems. XVOIGT allows the interactive fitting of Voigt pro- 
files to atomic absorption lines in astronomical spectra. We 
fitted simultaneously all identified lines in each system in or- 
der to determine the column densities and the Doppler widths 
of the systems, assuming that each system is formed by a sin- 
gle absorber. The large number of lines fitted in each sys- 
tem allowed us to fit even the saturated lines, using the dif- 
ferences in their saturated Voigt profiles. We find the —50 
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FIG. 3. — (a) Velocity spectra for the two normalized O VI lines. Air- 
glow geocoronal emission lines and interstellar Galactic absorption lines are 
marked with ©. The dashed line marks the rest frame zero velocity and the 
dotted-dashed line marks the recession velocity of NGC 4051, i.e., absorption 
around this velocity is suspected to be Galactic, (b) Velocity spectrum com- 
parison between the O VI A 1032 line and the N V A 1242 line from Collinge 
et al. (2001). The individual N V A 1242 absorption systems are blended into 
one wide trough in the O VI A 1032 line, (c) Velocity spectrum comparison 
between the C III A977 line and the C IV A 1550 line from Collinge et al. 
(2001). The flux of the C IV A1550 spectrum was scaled by 0.5 for clarity of 
the plot. There is a good correlation between the absorption systems seen in 
both ions, though the C III A977 spectrum has poorer resolution. Not all 10 
absorption systems seen in C IV A 1550 can be identified in the C III A977 
absorption, though the C III profile is consistent with all 10 systems being 
present. 

kms -1 system to have a column density of 2.2+Jj x 10 17 
cm~ 2 and the —240 kms -1 system to have a column density 
of l-O^Qj x 10 16 cm~ 2 ; we note that these column density 
estimates could have significant systematic errors if there is 
unresolved velocity substructure in the lines. For both sys- 
tems we find a Doppler width of b ~ 35 ± 15 kms" 1 . Short- 
ward of the Lyman edge the flux drops to about 30% of the 
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Velocity spectra for the first eight lines in the H I Lyman series. 

dashed line marks zero velocity. Airglow geocoronal emission 
erstellar Galactic absorption lines are marked with ©. The dashed 
he rest frame zero velocity and the dotted-dashed line marks the 
:locity of NGC4051, i.e., absorption around this velocity is sus- 
■ Galactic. 



continuum value. This can be explained by the total column 
density found above (~ 2.3 x 10 17 cm~ 2 ) which can produce 
such Lyman continuum absorption. The Galactic Lyman edge 
appears at 909.7 A in the NGC4051 rest-frame spectrum (see 
Figure^, and shortward of it the spectrum is totally absorbed. 

There is probably yet another absorption system of H I at 
around —400 to —500 kms -1 . This system is seen only in 
Ly/3, Ly5, Lye, and possibly Ly£ (Figure^. This absorption 
is shallower than the other two systems and is seen only in 
3-4 lines; hence we cannot measure its low column density as 
for the other two systems. 

Two Lya absorption systems with redshift velocities of 
~ 110 and ~ 260 kms -1 were suggested by Collinge et al. 
(2001) to originate in low-ionization, high-velocity clouds 
within NGC4051 (see their Figure 9). These are also iden- 
tified in Lyj3 absorption and marginally identified in Lyy and 
Ly<5 (Figure|4J and|4j). 

In the HETGS spectrum of NGC4051, Collinge et al. 
(2001) found a high-velocity absorption system at —2340 
km s~ 1 . We also looked for this system in the FUSE spectrum. 
For the O VI A 1032 line this high-velocity absorption system 
falls in the Lyj3 absorption trough where it is blended with 
geocoronal emission. For the O VI A 1038 line this absorption 
system falls on the blue wing of the O VI A 1032 trough. Fig- 
ure |3jt shows that the blue wings of the two O VI troughs are 
similar, although the blue wing for O VI A 1032 is a bit lower 
than for O VI A 1038. However, this difference is not signif- 
icant enough to claim a detection since such a small discrep- 
ancy could be caused by calibration effects or by uncertain- 
ties in the fit to the O VI emission. We also do not detect the 
—4500 kms~' absorption system reported by van der Meer 
et al. (2003) from the Chandra LETGS spectrum, nor do we 



detect the ~ —6500 kms or ~ —16500 kms (depending 
on the interpretation of an absorption line seen at ~ 7.1 keV) 
absorption system reported by Pounds et al. (2003) from an 
XMM-Newton observation. 

5. DISCUSSION 

5.1. Far Ultraviolet Emission and Absorption 

In these first FU V spectra of NGC 405 1 , we detect emission 
from O VI, N III, He II, and C III. The O VI emission lines are 
fitted with two Gaussians with FWHMs of ~ 700 and ~ 4000 
km s~ 1 ; the overall line profiles are significantly broader than 
for the optical lines. This is in agreement with the majority 
of other NLSls which show broad wings (several thousands 
of km s ) on their high-ionization UV lines (e.g., Rodriguez - 
Pascual, Mas-Hesse, & Santos-Lleo 1997). We also find the 
Gaussian centroids to be blueshifted by — 200kms~ 1 relative 
to the optical lines. This is again in agreement with the find- 
ings for other AGN (e.g., Richards et al. 2002 and references 
therein). 

We detect absorption on the blue wings of the O VI and C III 
emission lines, as well as absorption in the H I Lyman series. 
In H I we identify two main absorption systems and possibly 
a third at outflow velocities of -50 ±30, -240 ±40, and ~ 
—450 km s _1 . These three absorption systems are also present 
in C III (Figure 0;). The three FUV systems correspond in 
outflow velocity to UV systems 7-8, 5, and 2-3 which were 
found by Collinge et al. (2001) in N V, Si IV, and C IV. Thus, 
it is plausible that each of the three systems detected in the 
FUV is a blend of several narrower systems which cannot be 
resolved by FUSE? 

8 We note that UV system 5 (and perhaps other systems) of Collinge et al. 
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FIG. 5. — The combined 92.2 ks FUSE spectrum of NGC 405 1 binned to ~ 0.07 A pix 1 showing the Lyman absorption lines of H I toward the Lyman limit. 
All absorption lines up to Ly^ appear to be saturated and are blueshifted by —50 kms~'. The Galactic Lyman limit absorption starts at 909.7 A. The dashed 
horizontal line shows the zero flux level. Airglow geocoronal emission lines and interstellar Galactic absorption lines are marked with ©. 



An interesting relation is present between the strengths of 
the absorption systems and their outflow velocities. The ab- 
sorption system with the largest outflow velocity (~ —450 
kms -1 ) seems to have the weakest absorption, and it is 
present only in the lines up to about Ly£, i.e., it has a rela- 
tively low column density. The absorption system with the 
smallest outflow velocity (~ —50 kms -1 ) has the strongest 
absorption with relatively high column density (2.2+['i x 10 17 
cm -2 ), and all the lines that can be detected within the limit 
of the instrumental resolution (up to Ly^) are saturated. The 
intermediate-velocity absorption system (at ~ —250 kms -1 ) 
has an intermediate absorption depth, detected only up to 
LyA, and it has a column density of 1 .0^' 5 x 10 16 cm~ 2 . This 
relation might be explained if the systems with higher veloc- 
ities are closer to the central radiation source and hence are 
in a more intense radiation field which yields less H I in these 
systems. This might indicate that the outflow is decelerating 
so that its velocity decreases as its distance from the ioniz- 
ing source increases (a decelerating outflow was also found 
recently in NGC 3783; Gabel et al. 2003). In the lowest ve- 
locity outflow system, which is highly saturated, we detect on 
the trough's wing from 100 to —100 kms~' the effect of a 
velocity dependent covering factor (see Figure [3^), i.e., the 
absorption profile is produced by a different covering factor 
for each outflow velocity and not by the internal Voigt pro- 
files of the lines (e.g., Arav, Korista, & de Kool 2002). 

We can use the multi-epoch observations to put some con- 
straints on the line-of-sight acceleration of the absorbers. 9 

(2001) is saturated and broad; hence it might well be a blend of several sys- 
tems which are intrinsically blended (rather than just blended due to instru- 
mental resolution). 

9 This might not be the physical acceleration of the absorbing material, 
which could be crossing our line of sight to the nucleus in a "standing" pat- 
tern. 



First we use the FUSE observations which are separated by 
354.4 days in the rest frame. An upper limit on the absorption 
profile difference between the two epochs is 30 kms -1 (this 
is the trough location measurement uncertainty; here we used 
the O VI, C III, and Lyy to LyO absorption lines). Thus, the 
upper limit on the acceleration of the absorber is 9.8 x 10 
kms~ 2 . Next we used the time between the STIS observation 
and the last FUSE observation which is 1086.9 days in the rest 
frame. From the observed consistency of 100 km s between 
the two (e.g., Figure |3p), we derive an upper limit on the ab- 
sorber's acceleration of 1.1 x 10~ 6 kms~ 2 . The upper limit 
we find is at the lower range of the UV-outflow deceleration 
discovered recently in NGC 3783 of - (1-2.5) x 10~ 6 kms- 2 
(Gabel et al. 2003). 

The O VI absorption we detect from NGC 405 1 is an exam- 
ple of the "smooth" absorption morphology described in Kriss 
(2002). This morphology includes objects where the O VI ab- 
sorption is so broad and blended that individual O VI com- 
ponents cannot be identified. Another NLS 1 with a "smooth" 
O VI absorption morphology is Ark 564 (Romano et al. 2002). 
Other NLS Is have "single" morphology (like Ton SI 80 and 
IZwl) or a "blend" morphology (like Mrk478; see Kriss 
2002 for details). Thus, NLS Is appear to be similar in their 
range of O VI absorption-trough morphologies to broad-line 
Seyfert Is. 

The deep absorption troughs imprinted on the O VI and 
C III broad emission lines indicate these absorbers reside at 
a distance from the AGN center larger than the corresponding 
broad line region (BLR) size. In NGC 4051 the BLR size has 
been measured, using optical Balmer lines, to be 3 ± 1 .5 light 
days ([7.8 ± 3.9] x 10 15 cm; Shemmer et al 2003, see also Pe- 
terson et al. 2000). The BLR size for the UV lines (higher 
ionization lines) is known to be smaller than that of the op- 
tical lines (lower ionization lines) by a factor of — 2-3 (e.g., 
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Peterson & Wandel 2000). This is attributed to the fact that 
the ionizing flux density is higher closer to the central source. 
Taking these considerations into account, we estimate a lower 
limit for the distance of the FUV absorber from the central 
continuum source to be ~ 2 x 10 15 cm. 

The amount of nuclear FUV/UV light scattered around the 
absorbing material and into our line of sight is practically zero 
since the absorption troughs of the low-ionization ions (H I, 
C III, C IV, and N V) reach saturation at zero flux. We find 
the higher ionization ion O VI to have a covering factor of 
~ 0.9-0.95. This might indicate that the higher ionization ab- 
sorbers have a smaller covering factor than the low-ionization 
absorbers. 

5.2. Ultraviolet — X-ray Connection 

Several studies have suggested a link between the UV and 
X-ray absorbers in AGN (e.g., Mathur, Elvis, & Wilkes 1995; 
Crenshaw et al. 1999). The data presented in this paper enable 
us to examine this connection in NGC 405 1 . 

In Figure[6]we compare the N V A 1238 absorption from the 
STIS spectrum, the O VI A 1038 absorption from the FUSE 
spectrum, and combined line profiles of He-like ions and H- 
like ions from the Chandra HETGS data (Figure 3 of Collinge 
et al. 2001). The resolution and S/N of the HETGS spectrum 
do not allow precise measurement of the X-ray absorption 
systems. The velocity shifts and FWHMs of the FUV and 
UV absorption systems are consistent with the low-velocity 
(—600 kms~') X-ray absorber which is seen in both the He- 
like and H-like lines. 10 Both the He-like and the H-like pro- 
files seem to extend to higher blueshift velocities than the 
FUV and UV absorbers; however, given the HETGS resolu- 
tion and the poor S/N of the X-ray spectrum, it is not possible 
to determine the significance of this effect. We do not identify 
in the FUV and UV absorption lines the high-velocity absorp- 
tion system (at —2340 kms -1 ) observed in the X-ray spec- 
trum by Collinge et al. (2001). In particular, we do not detect 
this system in O VI. The high-velocity absorption system in 
the X-ray spectrum was detected only in the H-like ions of Si, 
Mg, Ne and O; it is onlymarginally seen in the He-like ions of 
these elements (FigurejfJ). Thus, it is not implausible that for 
the lower ionization ion O VI we do not detect this absorption 
system. Further observations, probably with Chandra, are re- 
quired to investigate the putative —2340 km s _1 absorption 
system further. 

McHardy et al. (1995) and Nicastro et al. (1999) mod- 
eled ROSAT observations of NGC 405 1 and claimed signifi- 
cant ionized X-ray absorption variability. Our observations, in 
contrast, fail to find any significant FUV absorption variabil- 
ity (though the limited temporal sampling might have caused 
a coincidence of observing on times when conditions were 
similar). This initially appears surprising if the UV and X- 
ray absorbers are connected. One possible explanation is that 
the low -resolution ROSAT data were confused by complex X- 
ray spectral variability; changes in the power-law continuum 
and/or soft excess could have created a false impression of 
ionized absorption variability. Alternatively, this could be at- 
tributed to FUV line saturation. Since most of the FUV ab- 
sorption lines we detect appear saturated, changes in the ab- 

10 We note that the profile of the low-velocity absorber has its main trough 
between ~ — 500 and —900 km s~ 1 in X-ray H-like lines, whereas in X-ray 
He-like lines it spans a wider range from ~ —100 to —900 kms~' (see Fig- 
ures|^J; andloTlV The more highly ionized absorber may be more concentrated 
at higher outflow velocities, or alternatively emission lines may be filling in 
the low-velocity part of the H-like absorption trough. 




-5000 -4000 -3000 -2000 -1000 



Velocity [km s ] 

FIG. 6. — Line velocity profiles, (a) The N V A1238 absorption from the 
STIS spectrum, (b) The O VI A 1038 absorption from the FUSE spectrum, 
(c) and (d) are the combined line profiles of He-like and H-like ions, respec- 
tively, from the Chandra HETGS data binned to 200 kms~' (Figure 3 of 
Collinge et al. 2001). The dashed vertical lines show 0, —600, and —2340 
km s~ 1 velocities. The wide absorption trough around —2000 km s~ 1 in panel 
(b) is from O VI A 1032. Airglow geocoronal emission lines and interstellar 
Galactic absorption lines are marked with ©. 



sorber might lead to little change in the absorption lines. A 
possible exception is the C III absorption which does not seem 
to be heavily saturated; the lack of apparent C III absorption 
variability supports the idea that the absorber does not change 
substantially on time scales of months. 

5.3. The Spectral Energy Distribution of NGC 4051 

In Figure[7]we show a spectral energy distribution (SED) of 
NGC 405 1 from the infrared to the X-ray. All data were cor- 
rected for the Galactic absorption (see caption for details), but 
no intrinsic reddening correction was applied. In the radio re- 
gion (which is not shown in the figure) the object is radio quiet 
with vF v « 10~ 16 -10~ 15 ergcm~ 2 s . The data we present 
are not contemporaneous, but overall most of the data repre- 
sent an average over time of the object's flux (see the figure 
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TABLE 4 

Luminosity of NGC 405 1 in Several Bands 



Wavelength [A] 


Data used 


Luminosity [10 40 erg s '] 


4000-7500 


Optical 


61 


1300-3300 


WE 


37 


1200-1700 


STIS 


11 


950-1150 


FUSE 


3.8 


6.2-24.8 


HETGS 


16 


1.9-6.2 


HETGS 


13 



NOTE. — Luminosity is computed using a cosmology of Hq = 70 
kms~' Mpc~', CIm = 0.3, and £1a = 0.7, which results in a luminosity dis- 
tance of 9.8 x 10 6 pc to NGC 4051. 



caption for details). The SED resembles the one presented 
in Komossa & Fink (1997), although here the SED is more 
detailed and has larger wavelength coverage. 

The mean spectrum of NGC 405 1 from the FUV to the op- 
tical is represented by data from four instruments and is av- 
eraged over timescales of years. The data from all four in- 
struments are in good agreement (notable is the agreement of 
the STIS spectrum which is at a consistent flux level with the 
averaged IUE spectrum observed one-to-two decades before). 
The optical continuum of NGC 405 1 varies little (e.g., Shem- 
mer et al. 2003), the UV continuum is also consistent with 
little variation, and the current work shows that NGC 4051 
did not vary significantly in the FUV. Nevertheless, there are 
notable variations in the optical line fluxes (e.g., Shemmer 
et al. 2003). Since the optical-to-FUV continuum does not 
vary much, a plausible conclusion is that the line variations 
are driven by the rapid and large soft X-ray and extreme UV 
variations and are not driven by the near or far UV as sug- 
gested before for this object. 

The observed SED of NGC 405 1 can be fitted with a power 
law (f v °c v a ) in limited spectral ranges. We find a to be 
-1.46 ±0.01 in the optical range (5500-7500 A) and -2.03 ± 
0.02 in the UV range (1160-3000 A). The FUV cannot be 
fitted with a power law as it is contaminated by emission and 
absorption, though Figure[7]suggests the UV slope continues 
into the FUV. In the X-ray range we find from the HETGS 
spectrum in the range 2-6 keV that a — —0.84 ±0.04. In 
Table|4] we give the luminosity of NGC 405 1 in several bands 
assuming Ho = 70 kms -1 Mpc~', £2m = 0.3, and Q.a = 0.7. 

Overall the SED shows a gradual rise in the emitted en- 
ergy from the hard X-rays to the far-IR. Most of the observed 
energy from NGC 405 1 is seen in the far-IR, while the FUV 
emission is about an order of magnitude weaker. The UV and 
FUV spectra do not show indications for a big blue bump. 
Such an observed SED could be explained by dust along the 
line of sight which reddens the optical-to-FUV spectrum and 
re-emits energy in the far-IR. This would be similar to the 
observed SED of Ark 564 in which significant reddening was 
found (Crenshaw et al. 2002; Romano et al. 2003). However, 
data from EUVE and ROSAT strongly constrain the amount 
of cold gas along the line of sight that could be associated 
with the putative dust (e.g., McHardy et al. 1995; Uttley et al. 
2000). One solution to this difficulty could be to invoke an 
ionized absorber with embedded dust (e.g., Brandt, Fabian, & 
Pounds 1996; Rraemer et al. 2000), although even in this case 
there is little evidence at present for X-ray absorption edges 
from the dust grains themselves (e.g., Komossa & Fink 1997). 
Also, the relatively high flux in the far-IR could have a signifi- 
cant contribution from the host galaxy (e.g., Ward et al. 1987; 
Elvis et al. 1994). 



Further evidence against a dust-reddening interpretation 
comes from consideration of the He II A1640/A4686 ratio. 
We estimate the observed ratio to be a surprisingly high 
19.2 ±6.3 from a STIS measurement of the He II A 1640 
line [(2.5 ±0.3) x 10~ 13 erg cm- 2 s" 1 ; Collinge et al. 2001] 
and a ground-based measurement of the He II A4686 line 
[(1.3±0.4) x lO^ergcm^s- 1 ; Peterson etal. 2000]. Such 
a value is not supportive of a reddening interpretation which 
predicts the ratio to be < 9 (compare with Crenshaw et al. 
2002). If the continuum is reddened then the absorption at 
1640 A is ~ 10 times larger than at 4686 A (by comparing 
the observed local-continuum flux ratio in Figure0to what is 
seen in a normal AGN). Such absorption would be expected 
to produce a He II A1640/A4686 ratio of ~ 1. While our 
measured line ratio should be treated with caution since vari- 
ability can add a factor of ~ 2 to the uncertainty, it still seems 
significantly higher than the expected ratio from reddening. 
More accurate and simultaneous measurements of this and 
other line ratios are needed to confirm this result. 

An alternative to models with significant reddening is that, 
perhaps due to its low black hole mass (Shemmer et al. 2003), 
NGC 405 1 may have a high-temperature big blue bump that is 
not apparent even in the FUV (compare our Figure0with Fig- 
ure 3 of Puchnarewicz et al. 1995). Discrimination between 
significant reddening and a high-temperature big blue bump 
is challenging, and we leave this task to a future study. For 
example, it will be necessary to determine unambiguously if 
the observed excess soft X-ray emission is due to a continuum 
soft X-ray excess (e.g., Collinge et al. 2001) or a relativisti- 
cally broadened O VIII recombination spectrum (Ogle et al. 
2004). 

6. SUMMARY 

Using three FUSE observations, we have presented the first 
measurements of the FUV spectral and variability properties 
of the intensively studied NLS 1 NGC 405 1 . Our main results 
are the following: 

1 . NGC 405 1 shows no significant FUV variability in ei- 
ther its continuum or line properties during three epochs 
spanning about one year. 

2. We detect FUV emission lines from O VI, N III, C III, 
and He II as well as blueshifted absorption lines from 
O VI, C III, and H I. 

3. The FUV absorption is generally coincident in veloc- 
ity with the UV absorption seen by STIS, although the 
lower FUV spectral resolution limits our ability to re- 
solve all of the absorption complexity seen in the UV. 
There is no evidence for new velocity components seen 
only in higher ionization FUV lines. 

4. The FUV absorption appears coincident in velocity 
with only the lowest velocity (—600 kms -1 ) X-ray ab- 
sorption system known. The several claimed X-ray ab- 
sorption systems at higher velocities are not detected in 
the FUV, although this does not rule out the existence 
of these systems. 

5. The H I Lyman series absorption lines are detected up 
to the Lyman edge, and two main velocity systems are 
seen in these lines. Column density estimates for these 
two systems indicate that the H I column density is 
lower for systems with larger outflow velocities. A pos- 
sible third velocity system is also seen. 
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FIG. 7. — The SED of NGC4051. Optical data are the average spectrum from Shemmer et al. (2003). IUE data are mean archive spectra taken over ~ 10 yr. 
The FUSE data are the average spectrum from this work. EUVE data are from Uttley et al. (2000), and the error bar represents the variability range. STIS and 
HETGS data are from Collinge et al. (2001); the HETGS data represent an average of the high and low states. The LETGS data are reported by van der Meer 
et al. (2003) and were taken while NGC4051 was in a lower X-ray flux state. The spectra were heavily binned for clarity. The infrared data were taken from 
the NASA Extragalactic Database (NED) and were obtained by IRAS in the far-IR and ground-based observatories in the near-IR. The data were corrected for 
the Galactic absorption, even though it is low. We used E(B — V) = 0.013 mag and the extinction curve of Cardelli, Clayton, & Mathis (1989) for the UV to IR 
spectra, and N H = 1.3 x 10 20 cm" 2 (Elvis, Wilkes, & Lockman 1989) with the WABS model in XSPEC (Morrison & McCammon 1983) for the X-ray spectra. 
No intrinsic reddening correction was applied. 



6. Any line-of-sight acceleration of the FUV absorption is 
constrained to be < 1 mm s~ 2 . 

7. The deep absorption troughs imprinted on the O VI and 
C III broad emission lines indicate that the FUV absorb- 
ing material resides at least ~2x 10 15 cm from the cen- 
tral continuum source. 

8. Several low-ionization ions have troughs that reach 
saturation at zero flux, significantly constraining the 
amount of nuclear FUV/UV light scattered around the 
absorbing material. 
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